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ABSTRACT 

The  stability  of  the  toroidal  portion  of  a  torispherical  shell  under 
internal  pressure  is  considered  from  the  point  of  view  of  the  linear 
buckling  theory.  A  detailed  stress  analysis  of  the  prebuckled  shell  is 
made  employing  asymptotic  integration.  The  change  in  potential  energy  of 
the  shell  is  then  minimized  using  a  Rayleigh-Ritz  procedure  for  actual 
computation  of  the  critical  pressure.  Numerical  results  reveal  that 
elastic  buckling  may  occur  for  very  thin  shells  whose  material  has  a 
relatively  high  value  of  the  ratio  of  yield  stress  to  elastic  modulus. 
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SYMBOLS 


a,0 

z 

A,B 

RjjR, 

u,v,w 


coordinates  of  middle  surface  of  shell 
coordinate  normal  to  middle  surface  of  shell 
Lame’’  coefficients  associated  with  a,  |3 
principal  radii  of  curvature 
displacements  in  the  a,  {3,  z  directions 

additional  displacements  in  the  meridional,  circumferential, 
and  normal  directions  after  buckling 


ea' W 

CTa,ap,CTap 


strains  at  any  point  in  the  shell 
stresses  at  any  poinl  in  the  shell 


U 

<p,e 

b 

a 


r 


strain  energy  of  thin  shell 

middle  surface  strains  and  curvature  changes  in  a  thin  shell 

angle  of  rotation  of  normal  to  middle  surface  about  tangent 
to  lines  (3  =  constant,  a  =  constant 

coordinates  in  meridional,  circumferential  direction  on 
toroidal  middle  surface 

radius  of  cross  section  of  torus 

distance  from  center  line  to  center  of  toroidal  cross  section 

r  =  a  +  b  sin  cp,  horizontal  distance  from  center  line  to 
point  on  toroidal  miudle  surface 


X 

X 

h 

E 

v 


b/a 


X  = 


2tt 


rr  -  2<Pr 


shell  thickness 
Young's  modulus  of  elasticity 
Poisson's  ratio 


flexural  stiffness  of  torus 


:  IJ 


P 


p 


applied  pressure 
critical  pressure 


Per 

n  wave  number  defining  number  of  circumferential  buckles  in 

torus 

VNe0  middle  surface  stress  resultants  prior  to  buckling 

Ve0  middle  surface  stress  couples  prior  to  buckling 


Iv 


INTRODUCTION 


Torispherical  shells  are  frequently  employed  as  end  closures  for 
cylindrical  shells  both  in  missile  design  and  in  a  wide  variety  of 
industrial-type  pressure  vessels.  Such  shells  generally  consist  of 
a  shallow  spherical  cap  joined  to  a  toroidal  segment,  Joined  in  turn 
to  a  cylindrical  shell.  This  combined  shell  is  then  subjected  to  in¬ 
ternal  pressure.  It  was  known1  that  large,  compressive  hoop  stresses 
were  developed  in  the  torus,  and  that  for  very  thin  shells,  elastic 
buckling  was  a  distinct  possibility.  For  shells  whose  thickness  is 
heavy  enough  to  avoid  buckling,  but  which  may  still  be  rather  thin 
compared  to  the  toroidal  radius,  plastic  deformation  may  occur.  This 
has  been  examined  in  detail  by  Drucker  and  Shield.8’®  However,  the 
elastic  stability  of  torispherical  shells  was  not  considered  a  matter 
for  concern  until,  recently,  such  a  shell  was  actually  observed  to 
buckle  under  internal  pressure.  (See  Figure  1.)  The  problem  is  an 
interesting  one,  since  the  prebuckling  state  of  stress  in  such  a  shell 
is  not  a  simple  one.  The  membrane  state,  often  employed  to  estimate 
prebuckling  stresses,  would  actually  predict  a  state  of  tension  in  the 
torus.  Consequently  a  thorough  stress  analysis  of  the  shell  prior  to 
the  onset  of  instability  must  be  made.  This  stress  analysis  is  also 
useful  (and  necessary)  in  ascertaining  the  dividing  line  between  those 
configurations  which  are  likely  to  buckle  and  those  which  are  likely 
to  undergo  plastic  deformation. 

In  this  report,  prebuckling  stresses  in  a  shell  of  the  type  shown 
in  Figure  2,  subjected  to  internal  pressure,  were  determined  by  asymp¬ 
totic  integration  techniques,  with  special  attention  being  devoted  to 
the  particular  solution.  The  vertical  support  on  the  cylinder  was  as¬ 
sumed  to  be  sufficiently  far  removed  from  the  Junction  of  the  torus 
that  its  specific  form  had  little  influence  on  the  toroidal  stresses. 
These  results  were  then  incorporated  into  the  stability  equations  for 
the  toroidal  segment,  and  numerical  results  were  obtained  by  applying  * 
a  Rayleigh-Ritz  approach  to  an  appropriate  potential  energy  expression. 
Since  no  prior  theoretical  results  appear  to  be  available,  an  experi¬ 
mental  investigation  of  the  same  problem  was  undertaken  simultaneously. 
The  results  of  the  two  investigations  compare  favorably.  Finally,  a 
somewhat  simplified  stability  criterion  is  proposed  which  provides  re¬ 
sults  of  essentially  the  same  degree  of  accuracy,  at  a  considerable 
saving  in  effort. 


ANALYSIS  OF  THE  ELASTIC  STASH. ,Y  OF  THIN  TOR  I SPHER I  CAL  SHELLS 
UNOER  UNIFORM  INTERNAL  PRESSURE 

Theoretical  Background 

Equations  governing  both  the  equilibrium  and  the  stability  of  a  thin 
elastic  shell  may  be  obtained  in  the  following  manner ,4»B  Beginning  with 
an  appropriate  set  of  strain -displacement  relations  and  a  stress -strain 
law,  one  forms  thejstrain  energy  integral  associated  with  the  thin  shell 
under  consideration.-  For  a  homogeneous  isotropic  thin  shell  this  may  be 
expressed  as 

U  "  |  JJf[Va  +  apep  +  V<ib]  AB  da  dP  dz  U) 

Next,  the  principle  of  virtual  work  and  the  principle  of  stationary  poten¬ 
tial  energy  are  applied.  In  this  connection,  consider  a  shell  in  a  state 
of  equilibrium  characterized  by  displacements  Uq,  Vq,  w0.  ThiB  state  is 
said  to  be  one  of  "neutral"  equilibrium  if  there  exists  an  adjacent  equi¬ 
librium  state  differing  from  the  first  by  an  infinitesimal  variation  in 
the  quantities  characterizing  that  state.  We  therefore  consider  an  alter¬ 
nate  equilibrium  state  characterized  by  displacements 

u  ■  u0  +  nUj  v  =  v0  +  riVj  w  “  w0  +  ,nwl  (2) 

where  ri  is  an  infinitesimal,  independent  of  a,  p. 

The  strain  energy  associated  with  this  second  state  becomes 

u  =  u0  +  njx  +  n3ua.  (3) 

The  principle  of  stationary  potential  energy  implies  that  for  the  second 
state  to  be  one  of  equilibrium, 

eu  +  6v  =  o,  (U) 

where  6V  is  the  negative  work  done  by  external  loads  due  to  a  variation  in 
displacements.  Considering  only  variations  in  ux ,  vx,  wx,  we  haves 

6U1  +  6Vl  +  ^(QJ,  ♦  6Va)  -  0.  (5) 

But,  since  the  first  state  was  assumed  to  be  an  equilibrium  state,  the 
principle  of  virtual  displacements  implies  that 

BUj  +  6Vt  =  0.  (6) 

Hence  the  condition  that  the  first  equilibrium  state  be  unstable  is  that 

6Ua  +  6Va  =  0.  (7) 

To  summarize,  equation  6  yields  equations  governing  the  equilib¬ 
rium  of  a  state  described  by  displacements  uG,  vOJ  wQ,  while  equation  7 
yields  equations  governing  the  stability  of  this  state.  This  general 


statement  of  the  problem  is  valid  for  both  large  and  small  deflections. 
Thus,  depending  upon  the  generality  of  tne  strain-displacement  relations 
and  the  stress-strain  laws  assumed,  a  variety  of  theories  is  possible. 
Kempner6  delineates  several  of  these  in  considerable  detail  and  with 
attendant  clarity. 


Toroidal  Geometry  and  St ra i n-D i tp 1 acement  Relations 

We  consider  the  coordinate  system  shown  in  Figure  2  with  9,0  the 
middle  surface  coordinates  in  the  meridional  and  circumferential  direc¬ 
tions  and  z  normal  to  both,  positive  when  directed  inward.  The  princi¬ 
pal  radii  of  curvature  R, ,  R^  are  given  by  R,  *  b  and  R,  sin  cp  =  a  +  b 
sin  <p  -  r.  The  so-callea  Lame*  coefficients  ror  this  system  are  Ax  M  b 
and  Aa  ■  Rs  sin  cp  *  r.  With  this  information  it  is  a  simple  matter  to 
obtain  from  Reference  U  or  5  (with  obvious  changes  in  notation)  the 
following  strain -displacement  relations,  valid  for  small  strains  but 
large  rotations: 

t  “A  +  i  (ii)e2) 
cp  cp  2  y 

«e  =  £0  +  |  (V) 

«cpe=  Ve  "  “¥"0  ^ 

where 


*<p  -  4*,  -  z*p  *e  *  *e  -  zxe 


=  ?  -  zX 

<P0  9e  ®0 


and 


y  u(%  ■ w) 


^(v,0  +  u  cos  cp  -  w  sincp) 


?«,e  ' 


(8b) 


-  V  cos  cp) 


(8c) 


\  *  -  ^e),<p  *e  -  i(>),e  -  ^£_(u)0) 


y<P0  =  E(T),cp 

(8d) 

“fy  “  “(w>6  +  v  sin  <P) 

“e  =  -  ^(w»cp  +  u) 

(8e) 

In  obtaining  these  relations,  in  addition  to  the  underlying  assumptions 
of  thinness  of  the  shell  and  validity  of  Kirchhoff's  hypothesis,  we 
have  made  the  further  assumptions  that  components  of  strain  and  rota¬ 
tion  are  small  compared  to  unity,  and  in  addition,  strains  are  small 
compared  to  rotations,  i.e.,  t  -  0(f)  -  0  (<d*)  «  1.  Further,  we  have 
assumed  that  ujj,  the  component  of  rotation  in  the  plane  of  the  middle 
surface  of  the  shell  is  much  smaller  than  the  other  two  components, 

(ik)  and  u>0.  Finally,  it  is  possible  to  effect  a  considerable  saving  in 
algebra  with  a  minimal  expense  of  accuracy  by  neglecting  in  u^,  and  Uty 
the  contribution  of  the  u,v  terms.  This  may  be  justified  on  an  order 
of  magnitude  estimate  of  the  terms  involved,  or  by  examination  of  the 
effect  of  this  approximation  on  numerical  results.  The  latter  rein¬ 
forces  the  former.  This  approximation  results  in  what  is  sometimes  re¬ 
ferred  to  as  a  Donnell-type  theory.  In  what  follows,  then,  we  take. 


-  ^(w,e)  “>0  =  -  £(w,cp)  .  (8f) 


With  a  stress-strain  law  assumed  in  the  form 


CT<P0 


2(l+v)  #CP0 


(9) 


the  strain  energy  of  the  shell  may  be  written  out  in  detail,  as  well  as 
the  quantities  and  U#,  sometimes  referred  to  as  the  first  and  second 
variation  of  the  strain  energy.  Specifically, 


ul  “  JXf  (°cD0efo1  +  <7eo,0l  +  of08o  )brdcod0dz 


(10) 


o 


(11) 


X  % 


u«  =  I  JIT  HV  +  *0^  +  Xec£ 

+  2[<*p0etp11  +  a0o«0ll  +  otp0o«peii  ])brdtpd0dz 


where 


V 

-  v  *  H  >* 

\  m  V  *  . )' 

V 

-  f*e  *» 

0 

0  u  0 

0  0^0 

0 

0  0  0 

V, 

•  \  *  w 

*e,  ’  \  *  'X\ 

Stp0l 

•  v,  -  XX 

-x\ 

V 

■  v 

x  -K)a 

Hi 

■  -  w  <12> 

and  where  it  ia  to  be  understood  that  the  subscripted  strain  or  rotation 
ia  to  be  evaluated  in  terms  of  the  displacement  of  the  same  subscript. 
Finally,  a y,  eg  and  ctyg  are  related  to  the  appropriately  subscripted 
strains  according  to  equation 


Stabi 1 ity  Criterion 

Employing  linearized  strain-displacement  relations,  we  may  obtain 
from  equation  6  the  equilibrium  equations  governing  the  (rotationally 
symmetric)  prebuckling  state  of  stress  in  the  shell.  We  prefer  to 
solve  these  in  a  somewhat  different  formulation  however,  and  omit  their 
presentation  here.  Proceeding  directly  to  the  stability  relations  we 
observe  first  that  the  prebuckled  state  is  one  of  small  deflections, 
governed  by  the  classical  (linear)  theory.  As  Novozhilov4  points 
out,  then,  we  are  justified  in  simplifying  the  expression  for  U„  by 
omitting  the  <a0  terms  in  ,  c0  and  e^g  .  Upon  integrating  through 
the  thickness,  h,  we  may  write  1  1 


"•  •  2(^7  JT  I2®1’ *  V  *  *  tAi, 

*  SS  j  V  +  V  +  2”\\  * 

+  |  J'J  |  NcpQ(  U)g^  )a  +  Ngo(^)a  -  2NCpgo(ui;pi'J)0i  )|brdcpd0 


brdcpd0 

brdcpd© 


(13) 


-  7  •• 


where  ,  Ng  are  the  stress  resultants  of  the  prebuckled  shell,  and 
due  to  rotational  symmetry  prior  to  buckling,  Nmo  *  0. 

The  potential  energy  of  the  external  force  system,  Va,  is  equal  to 
the  product  of  the  external  load  and  the  increase  in  volume  enclosed  by 
the  shell.  An  approximate  expression  for  Vg  for  the  toroidal  shell 
segment  under  uniform  internal  pressure  p  is8 


■  f  JT  j  "  w)  +  jK  v>0  +  u  C0SCP  -  w  3in  <p) 


vasin{P  ua 

— = - r  " 


vw 


>e 


<P 


brd<pd6 


(Hi) 


Stability  is  governed  by  the  relation 


6(Ua  ♦  V.)  -  0. 


(15) 


It  is  clear  that  solution  of  the  differential  equations  governing  the 
stability  of  the  toroidal  shell  segment  would  be  a  formidable  task  even 
if  the  prebuckled  state  were  a  simple  one.  For  this  reason  we  employ 
a  Raylelgh-Ritz  approach  and  obtain  an  upper  bound  for  the  critical 
pressure.  In  this  connection,  guided  somewhat  by  experimental  observa¬ 
tion  of  the  buckled  pattern  in  the  toroidal  segment,  we  choose  the  fol¬ 
lowing  set  of  displacements 

ux  M  A  ra  cos  X  (ep— «pQ )  cos  (n0) 

vx  ■  B  r3  sin  X  (9-<P0)  sin  (n8) 

wx  ■  C  ra  sin  X  (<p-<p0)  cos  (n8)  (16) 

where  n  is  the  number  of  lobes  in  the  circumferential  direction  and 

X,  5,  C  are  undetermined  coefficients.  The  variations  in  ux,  vx,  wx  are 
reduced  to  variations  in  A,  B,  C .  Inserting  equation  16  into  equation 
15  and  requiring  that  the  resulting  equation  be  satisfied  for  arbitrary 
values  of  6A,  6B,  6C,  we  obtain  a  system  of  three  linear  homogeneous 
algebraic  equations  in  A,  B,  C,  whose  determinant  set  equal  to  zero 
yields  an  estimate  of  the  critical  pressure.  Omitting  details  of  inter¬ 
mediate  algebra  we  may  write 


-  8  - 


(17a) 


Per  =  h  /  n\ 
E  2(l-us)b  \D  ) 

where 


N  *  Uauaaaa33  -  all(aB3)8  -  aia8  a33  +  a1#a13aa3  -  aaaa138 
D  a  ^a^a^h^  +  aa38  btl  +  a188  b33  +  2a13aaab13  +  al3a  baa 

"  aia(an3bi3  +  aia^»3^  "  ka33(a*8blx  +  aubaa) 
-KiaP8b3S  (17b) 

where i 


-  u  x»U-)(^)  -  uyXk^f)  ♦  ♦  n^) 

+  Uvn8^^-  2\3U+v)^^  +  2uX8n3^^^  -  2ns(l+3v)^^ 

♦  2n8X(l-6u)^+  2ns(l-v)  X8^^  J 


•'ll 


1  3 


Ja» 


“as 


■'a  a 


. .  i(m\ 

2\bB  / 

l/22*\ 

"2\b4  ) 

.  n^— ^ 

.  _  _  ±fe«\  +  jc  nfli  +  h 

2\bB  )  2\b *  )  ^  2  2t 


(17c) 


where  the  quantities  (1*),  etc., are  definite  integrals  of  the  form 

rr/2  _ 

(1*)  =  r  sina  X(  cp-cp0  )dcp, 

and  are  listed  in  the  Appendix.  All  these  integrals  may  readily  be 
evaluated  explicitly  in  terms  of  trigonometric  functions.  Finally, 

n/2  _ 

Ka  a  ^  Nq  r3sin8  X(  W-cp^dcp 

n/2  /  _  _  \  3 

Kx  ■  ^  Ncp  r^2r  r,^  sin  \(w-c p0)  +  X  r3  cos  X(«p-cp0)j  dqy. 

°  °  (18) 


Prebuckled  State  of  Strese 
In  Ut  we  encounter  the  integrals 

brdW>  i  II  ^fJbrdW  (19) 

which  assess  the  contribution  of  the  prebuckled  state  of  stress  to  the 
strain  energy  of  the  buckled  state.  Since  the  shells  under  consideration 
are  very  thin,  the  asymptotic  method  of  solution  of  Reissner's  equations 
for  symmetric  deformation  of  shells  of  revolution  is  most  appropriate. 


-  10  - 


Clark7*8  has  presented  asymptotic  representations  i'or  both  the  homoge¬ 
neous  and  particular  solutions  for  the  toroidal  shell.  These  were 
employed  in  the  present  investigation. 

In  this  connection,  it  may  be  well  to  digress  for  a  moment  to  dis¬ 
cuss  the  particular  integral  for  toroidal  shells.  Galletley1  has  given 
ample  warning  that  the  frequently  employed  practice  of  using  the  membrane 
solution  as  a  particular  solution  to  the  nonhomogeneous  thin  shell  equa¬ 
tions  is  not  valid  in  the  case  of  the  torus.  His  results  of  a  numerical 
integration  of  the  complete  differential  equations  for  just  such  a  tori- 
spherical-cylindrical  shell  show  a  markedly  different  stress  pattern 
when  compared  to  a  solution  using  the  membrane  state  as  a  particular 
solution.  However,  when  the  first  two  terms  of  an  asymptotic  particular 
solution  developed  by  Clark8  are  used  in  conjunction  with  the  homogeneous 
asymptotic  solution,  the  results  agree  remarkably  well  with  Galletley1 s. 
(See  Figure  3  •)  This  agreement  is  obtained  at  a  value  of  the  asymptotic 
parameter  (p  ■  15)  which  is  considerably  smaller  than  the  values  assumed 
by  the  shells  whose  stability  we  wish  to  consider,  and  thus  assures 
validity  of  application  of  the  asymptotic  method. 


Very  briefly,  the  prebuckling  stresses  may  be  written  in  terms  of 
the  two  basic  functions  J3  and  V  according  to  the  relations : 

N  =  ^-(¥  cos  cp  +  fl  sin  cp) 

<Po  rm 

"So  -fiX,  *^rbpH>  (20) 

where 

P  “  (l+Xsincp)-i/aQ  (A0hlr  -  Bphji  +  C0h8r  -  Dphai  +  |W3(\/D )i/». 

*  [F0Tr  -  G0Ti]  +(jfi»)-l(VU)l/aCG6(J"a  "  G(«p)D} 

¥  «  (1+Xsint pWa  Q{B0hir  ♦  A0hli  ♦  D0h2r  +  C0hHi+(kra/3( X/D )i/a* 

•  [F0Ti  +  G0Tr]) -  (pQJ)"1  (X/D)iA[F0W3  -  F(cp)3)  (21) 

and 


-  11  - 


(|®)i/e(c#,q))-i/a. 


n  "  H?  rV  "  "  Eh*  J  rbpVd<P^  F(<P}  ‘^h3'(i^sinip,))l/a(rV)coa  »* 

Q(tp)  7  I  (v  +  ^r2) (rV)  ‘ 2bSp  8ln  * cos  * 

-  bpr  cos  cp F0  b  F(tp0)  G0  =  G(co0),  (21a) 

Pjj  and  py  are  the  horizontal  and  vertical  components  of  applied  pressure. 
In  the  above,  hjr,  hji,  h8r,  hai  are  the  real  and  imaginary  components 
of  the  modified  Hankel  functions  of  order  one -third,  of  the  first  and 
second  kind  respectively.  Their  argument  is  understood  to  be  (iy). 

Ti  are  the  real  and  imaginary  parts  of  a  special  function  introduced 
by  Clark,  and  satisfying 


T"  -  iy  T(y)  =  1. 


(22) 


A0,  B0,  C  ,  D0  are  constants  to  be  determined  from  the  boundary  condi¬ 
tions  of  the  specific  problem  considered.  In  the  present  analysis,  a 
shallow  spherical  segment  was  assumed  joined  to  the  torus  at  cp  ■  cp0, 

and  a  cylinder  at  <p  *  The  vertical  support  on  the  cylinder  was 

assumed  to  be  sufficiently  far  removed  from  the  junction  of  the  torus 
that  its  specific  form  had  little  influence  on  the  form  of  the  toroidal 
stresses.  rfy  ,  Nqo,  and  Mqq  were  evaluated  along  the  length  of  the 

torus.  A  typical  set  of  results  is  shown  in  Figure  The  integration 
indicated  in  equation  18  was  carried  out  numerically.  When  this  is  done 
for  a  specific  choice  of  geometric  parameters,  pcr/E  remains  a  function 
of  n.  The  minimum  value  of  critical  pressure  was  found  by  calculating 
pcr/E  over  a  range  of  n,  subject  to  the  condition  that  n  be  an  integer. 


Simplified  Formula  for  Critical  Pressure 

The  procedure  outlined  above  was  carried  out  initially  using 
equation  17.  It  soon  became  apparent  that  certain  terms  could  be 
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omitted  without  significantly  affecting  the  numerical  results,  and  at 
the  same  time  affording  considerable  simplification  of  the  relation  for 
critical  pressure.  The  validity  of  the  approximation  rests  on  the  fact 
that  the  value  of  n  found  to  minimize  pcr/E  was  always  rather  a  large 

number  compared  to  unity.  This  result  is  consistent  with  experimental 
observations  of  the  buckling  pattern.  (See  Figure  1.) 

The  simplified  relation  for  critical  pressure  may  be  written: 


Per  h  N  ,  _ .  x 

T  "  2(i-?)bT  (n  >:>  1}  (23a) 

where 

5  “  1*  a33  -  an  5»3 

D  “  b33(^12  -  U  alt  a„)  (23b) 

M1 f)<¥) 

-  /l0*\ 

*88  “  V  b3  / 

-m) 

tt/2 

Ka  =  J*  N0  r3  sin3  X(cp-<p0)d<p 
'Po  ° 

K,  -  f/2  r(2rr>tD  sin  X(<p<p_)  +  X  r3(cosX(tp-<p0))fdcp. 
V  °  (23c) 
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Equation  23  is  to  be  preferred  over  equation  17  for  computational  pur¬ 
poses  when  n  is  large.  In  all  cases  considered  in  this  investigation, 
the  value  of  n  minimizing  pcr  was  greater  than  UO  and  in  most  cases  was 
greater  than  60. 


Numerical  Results  -  Comoarison  With  Experiment  -  Discussion 

Computations  of  the  critical  pressure  have  been  carried  out  thus 
far  for  a  limited  number  of  parameters.  A  typical  set  of  numerical 
results  is  shown  in  Figure  The  parameters  in  these  curves  correspond 
to  those  involved  in  an  experimental  program  concerned  with  the  same 
problem.  In  each  numerical  evaluation  of  the  critical  pressure,  a 
stress  analysis  of  the  unbuckled  shell  was  made,  the  parameters  Kj  and 
Kg  determined,  and  then  the  stability  criterion  evaluated.  A 

In  the  experimental  program,  scale  models  representative  of  those 
used  in  missile  applications  were  tested.  These  models  were  made  of 
poly -vinyl  chloride,  a  material  chosen  among  other  reasons  for  its  rel¬ 
atively  high  ratio  of  yield  strength  to  elastic  modulus.  Some  of  the 
results  of  these  tests  are  also  shown  in  Figure  5.  In  the  case  of  the 
thickest  shells  tested,  (h/b  >  .007)  the  disagreement  between  theory 
and  experiment  increased  significantly.  The  prebuckling  stress  analysis 
for  these  shells  revealed  that  at  pressure  levels  below  the  predicted 
buckling  pressure,  the  difference  in  principal  stresses  in  the  shell 
near  the  junction  of  torus  and  spherical  cap  had  exceeded  the  yield 
stress  of  the  material.  Since  our  analysis  has  assumed  elastic  behavior 
throughout,  it  would  not  be  expected  to  apply  to  the  experimental  mate¬ 
rial  in  this  range  of  the  parameters. 

The  following  experimental  result  may  also  be  of  some  interest.  An 
aluminum  torispherical  bulkhead  was  tested  under  internal  pressure.  The 
parameters  involved  were 


<P0  =35°  a  =  3U.U3"  b  =  18.07”  haverage  =  .081"  E  =  107  v  =  .3. 

Dimples  began  to  appear  in  the  toroidal  section  at  about  p  =  25  psi. 

The  pressure  level  was  increased  to  LtO  psi,  and  when  the  load  was  re¬ 
leased  the  dimples  remained.  The  theoretical  buckling  pressure  for  a 
shell  of  this  material  (  see  Figure  5)  is  bit  psi.  however,  the  prebuckling 
state  of  stress  reveals  that  in  the  vicinity  of  cp  =  U5°,  the  value  of 

Ocp-OQ 

— - —  on  the  inside  and  outside  surfaces  reaches  a  value  of  1300  and 
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1015,  respectively.  Thus,  for  aluminum  with  a  yield  stress  of  35,000 
psi  (as  was  the  case  in  the  test  described)  a  value  of  p  =  25  psi  would 
produce  stresses  very  close  to  yield. 

The  plastic  models  of  this  configuration,  however,  appeared  to 
buckle  elastically  at  3.2  psi.  (Theory  predicts  2.9  p3i.)  This  corre¬ 
sponds  to  an  elastic  buckling  pressure  of  70  psi  for  aluminum  with 
E  *  107 .  The  stress  level  induced  near  U5°  in  the  torus  in  the  plastic 
models  for  p  =  2.9  psi  is  in  the  vicinity  of  3700  psi,  which  is  below 
the  yield  stress  of  the  plastic  material  used. 

It  is  clear,  then,  that  the  phenomenon  of  elastic  buckling  of 
torispherical  shells  under  internal  pressure  occurs  only  for  very  thin 
shells  whose  material  has  a  relatively  high  value  of  the  ratio  of  yield 
stress  to  elastic  modulus.  With  the  increasing  role  being  played  by 
such  materials  in  space  structures,  it  is  believed  that  the  analysis 
described  in  this  report  will  assist  the  designer  toward  his  objective 
of  increased  efficiency.  Finally,  it  should  be  pointed  out  that  the 
analysis  may  be  applied  to  toriconical  shells  simply  by  replacing  the 
spherical  cap  by  a  cone  in  determining  the  prebuckled  state  of  stress. 
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Figure  I.  EXPERIMENTAL  SETUP 
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APPENDIX 


The  following  definite  integrals  occur  in  the  expression  for  criti¬ 
cal  pressure  (equation  17  or  23).  In  each  instance  the  limits  of  in¬ 
tegration  are  from  tpQ  to  tt/2.  The  variable  r  is  given  by  r  ■  a  +  b  sin  tp* 
All  integrals  may  be  evaluated  explicitly  in  terms  of  trigonometric 
functions . 

1*  ■  J  r  sin3  X(  tp-tp0  )dtp 
2#  a  J  r  sin  tp  cos3  tp  sin3  X(tp-tp0)dtp 
3*  *  J  r  cos3  tp  sin3  X(tp-tp0)dtp 

U*  ■  J  r  cos*  tp  sin3  X(tp-tp0)dtp 

5*  “  J  r  cos4  tp  sin3  X(tp-tp0)dtp 

6*  -  J*  r3  sin  tp  sin3  X(tp-tp0)dtp 

7*  ■  J*  ra  sin  tp  cos3  tp  sin3  X(tp-tp0)dcp 

8*  ■  j*  r3  cos  tp  sin  X(tp-tp0)cos  X(tp-tp0)dtp 
9*  ■  J*  ra  cos3  tp  sin  X(tp-cp0)cos  X(tp-tp0)dcp 
10*  “Jr3  sin3  X(te-tp0)dtp 
11*  “Jr3  cos3  X(tp-tp0)dtp 

12*  -  J*  r3  cos3  tp  cos3  X(tp-tp0)dtp 
13*  “Jr3  cos3  tp  sin3  X(cp-cp0)dtp 
ll»*  ■  (*  r3  sin3  tp  sin3  X(tp-cflQ)dtp 
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15* 

■  I  r“ 

cos3  ffl  sin3  \({p-5p0)d<P 

16* 

■  J r3 

cos  cp  sin3  X(<p-cp0  )dcp 

17* 

.  Jr3 

cos  <p  sin  X(<p-<P0)  cos 

X(cp-cp0)dcp 

18* 

,  jr3 

sin  9  cos  cp  sin  X(cp-cpc 

,)  cos  X(«p-tP0)dcp 

19* 

‘  J1* 

sin  <p  sin3  X(<p-tpo)d<p 

20* 

'  Jr‘ 

sin3  X(<p-cp0)dcp 

21* 

•  f  r4 

cos3  X(tp-cp0)d«p 

22* 

"  J  r4 

1  . 

sintp  sin3  X(cp-<p0)dcp 

23* 

“  f  r4 

sin  cp  cos3  X(cp-<p0  )dep 

2ii* 

■  J  r‘ 

cos3  <p  cos3  X(«p-cp0  )dcp 

25* 

cos  «p  sin  X(cp-^0)  cos 

X(cp-cp0)d«p 

26* 

.  Jr. 

sin3  X(  cp-cp0  )dcp 

27* 

*  J  r“ 

cos3  X(<p-w0)dcp 
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